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Abstract —  A  novel  adaptive  control  law  for  nonlinear  Hamil¬ 
tonian  Multi-Input  Multi-Output  (MIMO)  systems  with  uncer¬ 
tain  parameters  in  the  actuator  modeling  as  well  as  the  inertia 
and/or  the  Coriolis  and  centrifugal  terms  is  developed.  The 
physical  properties  of  the  Hamiltonian  systems  are  effectively 
used  in  the  control  design  and  the  stability  analysis.  The 
number  of  the  parameter  estimates  is  significantly  lowered 
as  compared  to  the  conventional  adaptive  control  methods.  A 
smooth  projection  algorithm  is  applied  to  keep  the  parameter 
estimates  inside  a  singularity-free  region.  The  developed  control 
scheme  is  applied  for  attitude  control  of  a  spacecraft  with  both 
the  inertia  and  the  actuator  uncertainties. 


I.  Introduction 

We  consider  a  multi-input  multi-output  (MIMO)  nonlinear 
Hamiltonian  system  represented  by  the  second-order  differ¬ 
ential  equation 

//(q)q  +  C(q,q)q  +  g(q)  =F  (la) 

F  =  D(q)u  (lb) 


where  q  £  R"  is  the  generalized  coordinates  vector,  H  £  R”x” 
is  the  (symmetric  positive  definite)  inertia  matrix,  Cq  is 
a  nonlinear  vector  of  Coriolis  and  centripetal  forces,  and 
g  £  R"  is  the  gravity  vector.  F  £  R”  is  the  generalized  force 
and  is  generated  by  a  control  input  vector  u  £  R”'  and  the 
actuator  matrix  D  £  R"xm.  For  full  tracking  control,  it  is 
generally  required  that  n  <  m  and  D  has  full  row  rank. 
Expressing  dynamics  of  systems  in  the  form  of  Eq.  (1), 
rather  than  the  state-space  form,  has  several  advantages. 
Equation  (1)  can  be  easily  derived  by  applying  Lagrange’s 
equation,  and  its  form  is  so  general  that  it  can  represent 
various  kinds  of  dynamic  systems,  such  as  a  multilink  robot 
manipulator  [1],  [2]  and  a  spacecraft  [3]— [5] ,  etc.  In  addition, 
there  is  a  physical  property  that  the  matrix  H  —  2C  is  skew- 
symmetric.  This  property  is  extremely  useful  in  designing 
advanced  control  schemes. 

Suppose  now  that  the  system  matrices  have  uncertainties 
in  their  parameters  and  can  be  expressed  as 


H(q,®s)=H"(q)+HA(q,®s)  (2a) 

C(q,q,0s)  =C”(q,q)+CA(q,q,0.v)  (2b) 

g(q,0.s)  =  g"(q)  +  gA(q,0s)  (2c) 

L»(q,0fl)=£»"(q)+£)A(q,0fl)  (2d) 


where  the  matrices  with  a  superscript  of  ‘n’  are  with  their 
known  nominal  values,  0fl  £  Rp  is  a  vector  of  unknown 
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bounded  constant  uncertainties  in  the  actuator  matrix  D, 
and  0V  £  W1  is  in  the  other  system  matrices/vector,  H,C, 
and  g.  We  also  assume  that,  by  proper  definition  of  the 
unknown  parameters  ®s  and  0a,  the  uncertain  matrices 
HA  (q,  ©i ) ,  CA  (q,  q,  ©5) ,  gA  (q,  0S)  and  DA(q.0„)  are  linearly 
dependent  on  0S,  and  ©a,  respectively. 

Adaptive  control  for  special  cases  where  the  actuator 
modeling  does  not  have  uncertainties,  that  is  Z)A  =  0,  has 
been  intensively  studied  in  the  literature  (see  for  instance 
Ref.  [1]).  However,  adaptive  control  for  more  general  cases 
with  Da  7/  0  does  not  seem  to  have  received  much  attention 
in  the  literature,  even  though  this  uncertainty  may  result  in 
significant  degeneration  of  controller  performance.  Ge  [6] 
has  derived  an  adaptive  control  law  for  multilink  manipulator 
systems  with  uncertainties  in  the  control  input  term,  but 
the  uncertainty  must  be  in  the  input  scalings,  and  thus  the 
uncertainty  matrix  must  be  diagonal  when  represented  in 
multiplicative  form.  (Or  it  can  be  said  that  DA  =  AD"  where 
A  is  a  diagonal  matrix.)  Chang  [7  ]  has  provided  an  adaptive, 
robust  tracking  control  algorithm  for  nonlinear  MIMO  sys¬ 
tems  which  is  based  on  the  “smooth  projection  algorithm,” 
which  has  also  been  used  in  [8]  and  [9]  for  adaptive  control 
of  SISO  systems.  More  recently,  one  of  the  authors  [10],  [11] 
has  also  provided  an  adaptive  control  scheme  based  on  the 
smooth  projection  algorithm  which  is  applied  to  spacecraft 
attitude  tracking  with  uncertain  misalignments/inertia  of  the 
actuator  flywheels.  However,  these  previous  results  [7],  [10], 
[11]  are  based  on  purely  mathematical  approaches  and  do 
not  exploit  the  useful  physical  properties  of  the  Hamiltonian 
systems.  More  significantly,  they  considered  MIMO  systems 
represented  by  the  differential  equations  without  any  terms 
in  front  of  the  highest  derivative  of  the  state  variable  vector, 
like  the  state-space  form.  Therefore,  in  their  methods,  Eq.  (1) 
would  need  to  be  converted  as 

q  =  -HlCq-H-1g  +  HlDu,  (3) 

in  which  the  uncertain  matrices  are  multiplied  to  each  other. 
Therefore,  in  order  to  design  adaptive  laws  based  on  the 
linear  dependency  of  uncertainties,  their  methods  need  a 
“over-parameterization”,  which  means  they  need  to  estimate 
the  combinations  of  the  elements  in  0a  and  0S,  and  thus  the 
number  of  parameter  estimates  would  significantly  increase. 

In  the  present  paper,  an  adaptive  control  algorithm  for  the 
general  case  where  the  uncertain  matrices  in  Eq.  (2)  are  all 
nonzero  is  developed  using  the  smooth  projection  algorithm, 
which  keeps  the  parameter  estimates  inside  a  properly  de¬ 
fined  singular-free  convex  set.  The  proposed  adaptive  law 
exploits  the  physical  properties  of  the  Hamiltonion  system 
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and  has  a  more  compact  form  with  a  smaller  number  of 
parameter  estimates.  The  proposed  method  is  then  applied  for 
a  spacecraft  attitude  control  problem  with  inertia/actuator  un¬ 
certainties.  Finally,  numerical  examples  with  the  spacecraft 
are  provided  to  validate  the  proposed  law. 

II.  Adaptive  Control  Law 

The  first  part  of  the  derivation  of  the  adaptive  control 
law  follows  the  standard  design  procedure  for  Hamiltonian 
systems  in  [1],  [4].  Let  us  assume  the  desired  trajectory 
(Id (t) ,  qd (t)  and  q d(t)  to  be  bounded.  The  tracking  error 
vector  is  defined  as  q  =  q  —  qd  and  a  measure  of  tracking  s 
and  the  reference  velocity  q,  are  defined  as 

s  =  <j  +  Aq  =  q-qr,  (4) 


and 

qr  =  q d  -  Aq,  (5) 

where  the  matrix  —A  is  assumed  to  be  Hurwitz.  Let  0*  be 
the  parameter  estimate  vector  and  let  0*  =  0*  —  0*  be  a 
parameter  estimate  error  vector,  when  *  is  a  or  s. 

As  suggested  by  Slotine  et.  al.  [1]  from  a  physical  insight 
that  qrHq  is  the  system’s  kinetic  energy,  the  following 
Lyapunov  function  candidate  is  defined  as 

v(t)  =  i  [sths  +  ©Jr;1©*  +  ©Jr;1©,]  (6) 

where  r„,rs  are  positive  definite  weighting  matrices. 

Differentiating  V(t)  with  respect  to  time  and  using  the 
skew-symmetry  property  of  the  matrix  H  —  2C  to  replace 
the  term  \sT Hs  with  sTCs,  one  can  have  the  following 
expression: 

V  =  sT(Du-Hqr-Cqr-g)  +  ©Jr;I6a  +  ©Jr;I0s.  (7) 


When  the  actuator  matrix  D  is  exactly  known  (i.e.,  DA  =  0) 
and  has  full  row  rank  (i.e.,  rank(D)  =  n),  one  can  easily 
design  an  adaptive  control  law  using  the  methods  proposed 
in  the  previous  works  [1]. 

However,  since  D  is  assumed  to  contain  unknown  param¬ 
eters  as  well  as  H,C  and  g,  a  novel  control  law  is  proposed 
as  follows: 


(Dn  +  Da)u  =  ( (Hn  +  Ha  )  qr  +  (Cn  +  CA )  q 

+(g"  +  gA)-^/s) 


r 

(8) 


where  the  matrices/vector  with  a  ‘hat’  symbol  are  constructed 
using  the  parameter  estimates  0„  and  0,  instead  of  the 
(unknown)  actual  parameters.  Kd  is  a  gain  matrix  which  is 
a  positive  definite.  When  the  matrix  Dn  +DA  is  assumed  to 
have  full  row  rank,  the  (weighted)  minimum  norm  solution 
is  given  by 


u  =  (Dn  +DaY(  {Hn+HA)q,-  +  (C”  +  CA)qr 


+  (g"  +  gA)-^s 


where  (■)  '  denotes  the  (weighted)  pseudo-inverse  of  a  matrix 

[5],  [11]. 


The  control  law  (9)  leads  to 

V  =  sT  |  {Dn  +  DA )  -  DA)  |  u  -  (Hn  +  Ha )  qr  -  (Cn  +  CA ) qr 

-  (g" + gA)]  +  ©J  r; 1  4  +  ©J  r; 1 4 

=  Sr  [(Hn  +  Ha) qr  +  (Cn  +  CA)qr  +  (g"  +  gA)  -  Kd s 
-  5Au  -  (Hn  +  Ha) q,  -  (< Cn  +  CA)qr  -  (g"  +  gA ) 

+  ©Jr; '©a  + ©Jr; '4 

=  sr  f/Aq,-  +  CAq,-  +  gA  -  Kd s  -  DAu 

+  ©Jr;l4  +  ©Jr;'4 

=  —sTKds  +  sT(HAqr  +  CAq,-  +  gA)  -  sTDAu 

+  ©Jr; 1 4  +  ©Jr; '  4 

(10) 

Since  the  uncertainty  matrices/vector  are  assumed  to 
depend  linearly  on  0’s,  one  can  define  known  regressor 
matrices  (in  fact,  row  vectors)  Ys  =  Tj(q,q,q,-,q,-,s)  €  Rlx? 
and  Ya  =  Ta(q,q,q,.,q,-,s,u)  e  Rlxp  such  that 

sT  (HAqr  +  CAq,  +  gA )  =  Ys®s  (11) 

and 

-sTDAu  =  Ya@a  (12) 

Notice  that  the  definition  of  the  regressor  matrix  Ys  in 
Eq.  (11)  is  slightly  different  from  that  in  the  previous 
works  by  Slotine  [1],  [3].  By  including  s  in  the  regressor’s 
definition,  the  size  of  the  regressor  Ts  becomes  smaller  than 
that  of  Slotine’s  works  where  the  regressor  has  a  size  of 
n  x  q. 

The  time  derivative  of  V  then  becomes 

v  =  —sTKds + Tv0.v  +  L;4  +  ©Jr, 1 4  +  0 J r„ 1 4  (13) 


and  taking  the  adaptation  laws  of  the  parameter  estimates  to 
be 


©>■ 

II 

P 

(14) 

and 

4  =  -r aYj 

(15) 

then  yields 

V(t)  =  -sTKds  <  0. 

(16) 

Using  standard  arguments  in  [1],  [11]  which  use  Barbalat’s 
lemma,  one  can  easily  show  that  V  — >  0  and  thus  s  — >  0  as 
t  — >  °°.  This  also  implies  the  tracking  error  q  — >  0  as  well. 

A.  Smooth  Projection  Algorithm 

We  previously  used  an  assumption  that  the  matrix  Dn  +DA 
has  full  row  rank  in  deriving  the  adaptive  control  laws, 
Eqs.(9)  ,(14)  and  (15).  In  general,  for  the  full  tracking 
control,  the  nominal  matrix  Dn  in  general  has  full  row  rank. 
However,  a  drift  of  the  parameter  estimates  0„,  governed  by 
an  update  law  (15),  can  result  in  D"+DA  losing  rank.  We 
will  refer  to  this  situation  as  a  “singularity”  of  the  steering 
law  due  to  the  adaptation.  This  singularity  hinders  the  use 
of  the  derived  control  laws,  so  they  need  to  be  modified. 
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If  the  nominal  matrix  D"  has  full  row  rank,  and  the 
true  value  of  the  parameter  uncertainty  &a  is  bounded  by 
a  sufficiently  small  number,  and  the  parameter  estimate  0,, 
is  also  kept  small,  then  the  matrix  D"  +  DA  will  also  have  full 
row  rank.  To  this  end,  we  define  the  following  two  convex 
sets, 


we  choose  the  initial  parameter  guess  0„(O)  inside  the  set 
n(.)n .  For  instance,  we  may  take  0n(O)  =  0. 

It  is  also  worth  mentioning  that  the  proposed  adaptation 
law  (19)  has  the  additional  benefit  of  keeping  the  parameter 
estimates  from  “bursting”,  which  may  happen  when  the 
persistency  of  excitation  condition  does  not  hold  [14]. 


^®fl  =  {0aGKp|  ||0u||2</3}  (17) 

^  =  {0«er|  \\®a\\2<p  +  S}  (18) 

where  /3  >  0  and  <5  >  0  are  known  constants.  Notice  that 
C  £2©a .  We  make  the  following  three  assumptions. 

.  Assumption  1.  The  nominal  value  D"  has  full  row  rank 
of  n. 

•  Assumption  2.  The  actual  value  0„  belongs  to  the  set 

•  Assumption  3.  If  0a  £  £2©a,  then  D"  +  DA  is  non¬ 
singular. 

These  assumptions  allow  us  to  modify  the  adaptation  law 
(15)  by  using  the  “smooth  projection  algorithm”  as  follows1: 

4,=Proj(0„,<Ffl)  (19) 

where 

<&a  =  -r aYj  (20) 


and 


if  (i)  ||0a||2  <  p  ,or 

(ii)  1 1  ©«  1 1 2  >  P  and  <Fj0a  <  0, 


Froj  ( ©a .  <t>a ) 


flu  ■ 


|0„||2-j3)<I >l&a 


0, 


4||©fl||2 

if(iii)  ||®fl||2  >  P  and  ^Ta®a  >  0. 


(21) 

This  adaptation  law  is  identical  to  (15)  in  cases  (i)  and  (ii), 
and  switches  smoothly  to  a  new  expression  in  case  (iii). 
The  projection  operator  Proj (©a,<I>0)  is  locally  Lipschitz  in 
(0a,  <!>«),  thus  the  system  has  a  unique  solution  defined  for 
some  time  interval  [0,T),  T  >  0. 


Proposition  1:  Under  Assumptions  1,2,  and  3,  the  control 
law  Eq.  (9)  along  with  the  adaptation  laws  Eqs.  (14)  and 
(19)  yields 

V  <  -sTKds  <  0,  (22) 

and 

®a(t  =  o)  G  =>  &a(t)  e  Cl&a,  Vf  >  0.  (23) 

Proof:  The  proof  is  straightforward  and  therefore 

omitted  here.  It  is  similar  with  the  proof  in  the  author’s 
previous  works  [11],  [13],  ■ 

From  Proposition  1,  one  can  conclude  that,  using  the 
feedback  control  law  (9)  and  the  adaptation  laws  (14)  and 
(19),  q  — >  0  as  t  — >  °°  and  ( D "  +DA)  will  not  lose  rank,  if 

1  The  adaptation  law  is,  in  fact,  only  Lipschitz  continuous,  not  continu- 
ously  differentiable.  The  use  of  the  term  “smooth”  is  a  slight  misnomer  in 
this  context,  but  we  use  it  here  in  accordance  to  prior  usage  in  the  literature. 

It  should  be  noted  that  a  new  parameter  projection  operator  which  is 

has  been  recently  introduced  in  Ref.  [12]. 


III.  Application  to  Spacecraft  Attitude  Control 
A.  Equations  of  Motion 

In  this  section,  applying  the  proposed  adaptive  control 
scheme,  we  design  an  adaptive  attitude  tracking  control  law 
for  a  spacecraft.  A  cluster  of  Variable-Speed  Control  Moment 
Gyros  (VSCMGs)  with  N  flywheels  is  used  for  the  torque 
actuator.  While  a  conventional  Control  Moment  Gyro  (CMG) 
keeps  its  flywheel  spinning  at  a  constant  rate,  a  VSCMG  -  as 
its  name  implies  -  is  essentially  a  single-gimbal  CMG  with 
the  flywheel  allowed  to  have  variable  speed.  (See  Refs.  [10], 
[15]  for  more  details  and  applications  of  VSCMGs.) 


Fig.  1 .  Spacecraft  Body  with  the  i-th  VSCMG. 

Figure  1  shows  a  spacecraft  with  the  i-th  VSCMG,  where 
gi  is  (body-fixed)  gimbal  axis,  s,  is  spin  axis,  and  t,  =  g,  x 
s,  is  transverse  torque  axis.  The  equations  of  motion  of  a 
spacecraft  with  VSCMGs  are  complicated  as  shown  in  the 
aforementioned  references,  but  under  assumptions  which  are 
standard  in  the  literature  [11],  [13],  they  can  be  simplified 
as  follows: 

Jco—  [hx]m  +  gu  =  0,  (24) 

and 

h  =  J(o+AsIwCl,  (25) 

where 

Q  =  [Qc,ng,Qm\e  R3x2N,  (26) 

and  where  Qcmg{yM)  =  AtIw£ld ,  Qrw(Y)  =  ASIW,  and  the 
control  input  of  this  system  is 

u  =[ul,...,u2N)T  =  [fAT]TGR2N.  (27) 

J  is  the  total  moment  of  inertia  of  the  spacecraft  which  is 
assumed  to  be  constant2,  co  is  the  body  rate  vector  of  the 
spacecraft,  h  is  the  total  angular  momentum  of  the  spacecraft, 
y  =  [yi,...,yv]7’  e  and  Q.  =  [Qi , . . . , Q.N]T  e  Rw  are 
vectors  of  gimbal  angles  and  flywheel  spinning  speeds, 

2 In  fact,  J  is  a  function  of  y  but  the  dependence  is  weak  in  general. 
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respectively,  and  Iw  is  a  diagonal  matrix  with  the  inertias  of 
VSCMGs  flywheels.  The  skew-symmetric  matrix  [vx],  for 
vSR3,  represents  the  cross  product  operation.  The  matrices 
A*  G  R3xiv  have  as  columns  the  gimbal  (g,),  spin  (s,)  and 
transverse  (t,)  directional  unit  vectors  expressed  in  the  body- 
frame,  where  *  is  g,s  or  t.  These  matrices  depend  on  the 
gimbal  angles  as  follows 

As  =  Ago  (28) 

A.s  =  As0[cos7]6,+Ato[siny]'/  (29) 

A,  =  At0  [cos  yf-As0  [sin  yf  (30) 

where  the  A*o’s  denote  the  values  of  A*  at  7=  0.  The  symbol 
\d  denotes  the  diagonal  matrix  with  elements  the  components 
of  the  vector  x,  and  cos 7=  [cosyi,---  ,cosyv]r  and  sin 7  = 
[sin  7i ,  •  ■  ■  ,sinTv]7'. 

The  modified  Rodrigues  parameters  (MRPs)  [  1 6]— [  1 8]  are 
chosen  to  describe  the  attitude  kinematics  of  the  spacecraft.3 
The  kinematics  in  terms  of  the  MRPs  is  given  by 

a  =  G(a)co,  (31) 

where 

G(g)  =  ^I3  +  [c7x]  +  aCTr-[i(l  +  a7'ct)]l3)  (32) 

and  lr  is  the  r  x  r  identity  matrix. 

As  suggested  in  Refs.  [1]  and  [5],  we  combine  the  kinetic 
equation  (24)  and  the  kinematic  equations  (31)  into  one 
second-order  system  as  follows: 

H(o)o  +  C(g,o)g  =  D(g)u  (33) 

where 

H(g)=G~tJG~\  (34) 

C(g,g)  =  —G~TJG~lGG~1  —  G_r[(R(a)h7)x]G_1,  (35) 

D(g)  =  —G~tQ.  (36) 

R(g)  is  a  rotational  matrix  from  the  inertial  frame  to  the 
body  frame,  and  h7  is  the  total  angular  momentum  of  a 
spacecraft  expressed  in  the  inertial  frame  which  is  conserved 
to  be  constant  if  there  is  no  external  torque  applied  to  the 
spacecraft.  Therefore,  h  =  R(a)h7. 

Notice  that  the  equation  of  motion  (33)  has  the  form  of 
(1)  with  the  gravitational  term  g  =  0.  Moreover,  it  can  be 
easily  shown  that  the  matrix  H  —  2C  is  skew-symmetric  [1]. 

B.  Adaptive  Attitude  Tracking  Control 

Suppose  that  there  are  uncertainties  in  D  as  well  as  J.  We 
assume  that  the  exact  values  of  the  initial  axis  directions  of 
VSCMGs  actuator  at  7  =  0  are  unknown.  This  can  happen 
when  the  VSCMGs  are  installed  with  small  misalignments 
and/or  the  measure  of  gimbal  angles  has  constant  unknown 
bias.  In  addition,  h7  is  also  unknown  constant  not  only 
because  of  uncertain  J  but  also  because  of  uncertain  As. 

3We  hasten  to  point  out  that  the  use  of  the  MRPs  to  describe  the 
kinematics  is  done  without  loss  of  generality.  Any  other  suitable  kinematic 
description  could  have  been  used  with  the  conclusions  of  the  paper  remain¬ 
ing  essentially  the  same. 


For  most  cases  the  effect  of  axes  uncertainties  on  the 
overall  system  performance  is  not  significant.  However,  for 
the  case  of  flywheels  used  as  “mechanical  batteries”  in  an 
Integrated  Power  and  Attitude  Control  System  (IPACS)  [5], 
[19],  [20],  even  small  misalignments  of  the  flywheel  axes 
can  be  detrimental.  Flywheels  for  IPACS  applications  spin  at 
high  speeds  and  have  large  amounts  of  stored  kinetic  energy 
(and  hence  angular  momentum).  Precise  attitude  control 
requires  proper  momentum  management,  while  minimizing 
spurious  output  torques.  This  can  be  achieved  with  the  use 
(in  the  simplest  scenario)  of  at  least  four  flywheels,  whose 
angular  momenta  have  to  be  canceled  or  regulated  with 
high  precision.  If  the  exact  direction  of  the  axes  (hence 
the  direction  of  the  angular  momenta)  are  not  known  with 
sufficient  accuracy,  large  output  torque  errors  will  impact  the 
attitude  of  the  spacecraft. 

The  uncertain  parameters  in  H  and  C  can  be  defined  as 
follows: 

©,  =  [Ajll,Aj22,Aj33,Aji2,A/i3,A723,A/l1,A/l2,A/?3]:r  G  R9 

(37) 

where  A/’s  are  the  elements  of  JA=J  J"  and  J”  is  the 
nominal  value  of  the  actual  inertia  matrix  J.  Similarly,  Ah’s 
are  elements  of  h7A  =  h7  -  h7'1 .  The  uncertain  parameters  in 
D  are  defined  as 


0a  =  [@,ri  >  •  •  •  ,  ©iyv,  ®£l ,  •  -  .  ©£#] 1 


j)6Af 


(38) 


where 


©;,<  —  t;,0  -  t"0,  ®i, i  —  S;,0  -  S" 


0) 


1=1, 


,N.  (39) 


and  t,  o  and  S/.o  are  actual  value  of  t,  and  s ,•  at  7  =  0, 
respectively,  and  t"0  and  s"0  are  their  nominal  values.  The 
total  number  of  parameter  estimates  is  then  6/V  +  9.  If  one 
uses  the  methods  in  the  previous  works  [7],  [10],  [11], 
[13],  then  the  number  of  the  estimates  will  be  as  much  as 
6  x  (6/V +  3). 

Exact  mathematical  expressions  of  the  regressor  vectors  Ys 
defined  in  (11)  can  be  easily  obtained  using  symbolic  math 
packages,  and  can  be  constructed  from  the  measurements  of 
a,  <7,  and  the  desired  trajectories  <7j,  07/,  Oj.  The  regressor 
vector  Ya  defined  in  (12)  can  be  obtained  in  the  same  way, 
but  it  is  also  possible  to  derive  its  mathematical  expression 
by  manipulating  the  matrices  as  follows: 

Ya  =  sTG-T[(Mu)lh,---,(Mu)2NI3]  (40) 


where 


[cosy]'7 

[siny]'7 

'  iHnd  o 

—  [siny]'7 

[cosy]'7 

0  4  _ 

(41) 


and  (x),-  is  the  ith  element  of  a  vector  x.  Then  using  the 
developed  control  scheme  in  Sec.  II,  one  can  design  an 
adaptive  attitude  tracking  control  law  for  the  spacecraft. 


IV.  Numerical  Examples 

Numerical  examples  for  a  satellite  with  a  VSCMGs  cluster 
are  provided  in  this  section  to  test  the  proposed  adaptive 
control  algorithm.  A  standard  four-VSCMG  pyramid  con¬ 
figuration  ( N  =  4)  is  utilized  [20].  The  parameters  used  for 
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TABLE  I 

Simulation  Parameters  and  Gains 


Symbol 

Value 

Units 

<7(0) 

[0,0,0]' 

- 

ffl(0) 

[0,0, 0]r 

rad/sec 

ffl(o) 

[0,0, 0]r 

rad/ sec2 

7(0) 

[0,0,0,0]r 

rad 

7(0) 

[0,0,0,0]r 

rad/sec 

£2(0) 

104  x  [2.5,3.5,3.5,3.0]r 

rpm 

Av 

diag{2.0,  2.0,  2.0,  2.0} 

kgm2 

Kd 

103I3 

- 

A 

I3 

- 

A 

0.01 

- 

5 

0.01 

- 

the  simulations  are  shown  in  Table  I.  Notice  that  the  initial 
wheel  speeds  of  the  VSCMGs  are  set  to  25,000~35,000 
RPM,  which  are  an  order  of  magnitude  larger  than  the  speed 
of  conventional  CMGs,  since  the  flywheels  of  VSCMGs 
used  for  IPACS  in  general  need  to  spin  very  fast  so  that 
they  are  competitive  against  traditional  chemical  batteries. 
According  to  Ref.  [21J,  even  a  higher  speed  than  these  vales 
is  implementable,  at  least  in  a  laboratory. 

The  nominal  values  of  the  axis  directions  at  7=  [0, 0, 0, 0]r 
are 


'  0 

-1 

0 

1  ' 

11 

1 

0 

-1 

0 

7 

(42) 

0 

0 

0 

0 

-0.5774 

0 

0.5774 

0 

0 

-0.5774 

0 

0.5774 

•  (43) 

0.8165 

0.8165 

0.8165 

0.8165 

The  (unknown)  actual  axis  directions  at  7  =  0  used  in  the 
present  example  are  assumed  as 


-0.0072 

-0.9999 

-0.0072 

0.9999 

0.9999 

0.0143 

-0.9999 

0.0057 

5 

-0.0071 

0.0071 

-0.0142 

-0.0117 

(44) 

-0.5657 

-0.0024 

0.5831 

0.0061  ' 

0.0018 

-0.5715 

-0.0157 

0.5868 

5 

(45) 

0.8246 

0.8206 

0.8123 

0.8097 

which  are  obtained  by  rotating  each  nominal  gimbal  frame 
(g;,s,',t,')  with  1  degree  about  arbitrary  direction.  With  these 
values,  1 1 ©„  1 1 2  ~  0.0015.  The  nominal  value  of  the  spacecraft 
inertia  matrix  is 


15000  3000  -1000 


Jn  = 

3000 

6500 

2000 

-1000 

2000 

12000 

(unknown)  actual  inertia  matrix 

18000 

2400 

-1200 

J  = 

2400 

5200 

1600 

-1200 

1600 

14400 

kg  m2 


kg  nr 


(46) 


(47) 


reference  attitude  is  aligned  with  the  body  frame  which  is 
also  aligned  with  the  inertial  frame,  and  the  angular  velocity 
of  the  reference  attitude  is  chosen  as 

,  .  .  2m  .  2nt  .  2 nt\T 

m{t)  =0.04 sin— , sm  — , sin  —  J  rad/sec.  (48) 

First,  in  order  to  show  the  effect  of  the  misalignment  of  the 
axes  of  the  VSCMG  cluster  and  the  uncertain  inertia  matrix, 
a  simulation  without  adaptation  was  performed.  Figure  2 
shows  the  attitude  tracking  error  (expressed  with  ‘3-2-1’ 
Euler  angles)4  under  control  law  with  the  adaptation  gains 
rfl  and  rs  set  to  zero  matrices.  Since  the  flywheel  speeds 
are  very  fast,  there  is  large  attitude  tracking  error  without 
adaptation. 

Next,  another  simulation  was  run  with  adaptation  of  the 
actuator  uncertainty  0a  only.  The  adaptation  gains  are  set 
to  Fa  :  1 0  I6/V  and  rs  =  0,  and  the  resulting  attitude 

error  is  shown  in  Fig.  3.  There  is  significant  performance 
improvement  with  the  adaptation  of  0„  only,  but  tracking 
error  with  a  magnitude  of  about  0.2  degree  remains.  On 
the  other  hand.  Fig.  4  shows  the  attitude  tracking  error  with 
adaptation  of  0S  only.  The  adaptation  gains  are  Fa  =  0  and 
rs  =  diag(107  l6,105  I3).  In  fact,  the  control  law  in  this 
scenario  is  almost  identical  with  Slotine’s  method  [1],  The 
attitude  error  is  significantly  attenuated  using  the  adaptive 
controller,  but  there  are  again  residual  tracking  errors  with  a 
magnitude  of  about  0.1  degree. 

Finally,  a  simulation  is  performed  with  adaptation  of  both 
Ta  and  that  is  =  10  Ig/y  and  rs  =  diag(107  Ig ,  105  I3). 
Figure  5  shows  the  tracking  performance  is  improved  upon 
Slotine’s  control  law.  Figure  6  shows  the  time  history  of 
||0fl||2.  It  is  confirmed  that  ||0a||2  does  not  drift  more  than 
ft  +  8  =  0.02  owing  to  the  smooth  projection  algorithm.  As 
a  result,  the  steering  law  (9)  remains  well-defined. 

V.  Conclusions 

In  this  paper,  we  proposed  an  adaptive  tracking  control  law 
for  a  nonlinear  Hamiltonian  MIMO  dynamic  system.  The 
proposed  control  scheme  has  several  significant  improve¬ 
ments  over  the  previous  works  in  the  literature.  First,  the 
proposed  method  can  deal  with  uncertainties  in  the  actuator 
terms  which  Slotine’s  method  [  1  ]  does  not  deal  with.  The  re¬ 
gressor  matrix  for  adaptation  of  the  uncertain  inertia  also  has 
a  smaller  size  than  that  of  [1].  Second,  the  proposed  method 
exploits  the  physical  properties  of  the  Hamiltonian  systems 
and  so  the  designed  law  is  more  compact  than  those  in  [7], 
[11],  [13]  which  deal  with  actuator  uncertainties  but  are 
derived  based  on  purely  mathematical  approaches.  Finally, 
the  proposed  method  deals  with  the  actuator  uncertainties 
separately  from  the  uncertain  inertia/Coriolis/gravity  terms. 
Therefore,  it  does  not  need  over-parameterization  to  deal 
with  both  kinds  of  uncertainties  at  the  same  time,  while  [7], 
[11],  [13]  do. 

The  developed  adaptive  algorithm  is  shown  to  significantly 
improve  the  tracking  performance  in  the  application  to  the 


which  is  obtained  by  adding/subtracting  20%  of  the  nominal  4The  use  of  Euler  angles  in  the  figures  in  done  solely  for  the  convenience 
values.  The  reference  trajectory  is  chosen  so  that  the  initial  of  the  reader  who  may  not  be  familiar  with  the  MRPs. 
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spacecraft  attitude  control,  but  it  still  has  room  for  improve¬ 
ment.  For  instance,  while  only  3  parameters  are  generally 
needed  to  express  a  misalignment  of  axis  frame  for  one 
VSCMG,  a  total  of  6  parameters  are  used  in  this  method. 
Development  of  methods  to  reduce  the  number  of  estimated 
parameters  would  be  extremely  beneficial. 
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Fig.  2.  Tracking  Error  Without  Adaptation. 

Attitude  Error 
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Fig.  3.  Tracking  Error  With  Adaptation  of  ®a  Only. 
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Fig.  4.  Tracking  Error  With  Adaptation  of  &s  Only. 
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Fig.  5.  Tracking  Error  With  Adaptation  of  0fl  and  ©5. 
Square  of  Norm  of  Parameter  Estimates 
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Fig.  6.  Square  of  Norm  of  Parameter  Estimate  ©fl. 
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